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Abstract.  In  order  to  obtain  optimal  performance  from  a free  electron  laser  (FEL),  the  axial  velocity 
spread  on  the  electron  beam  must  be  small  as  it  propagates  through  the  wiggler  field.  Treated 
separately,  both  the  wiggler-induccd  betatron  motions  and  the  self-induced  spaee  charge  forees  tend 
to  increase  the  axial  velocity  spread  and  degrade  the  performance  of  the  FEL.  However,  it  has  been 
shown  analytically  [B.  Hafizi  and  C.W.  Roberson,  Phys.  Plasmas  3,  2156  (1996)]  that  when  both 
effects  are  treated  self-consistently  an  equilibrium  exists  wherein  the  space  charge  forces  exactly 
compensate  for  the  betatron  motion.  This  leads  to  the  surprising  result  that  for  a continuous  beam, 
increasing  the  beam  current  can  improve  the  beam  quality. 


INTRODUCTION 

The  quality  of  the  electron  beam  used  in  a free  electron  laser  (FEL)  can  be  quantified  in 
terms  of  the  scaled  thenual  velocity,  5,  given  by  [1] 

5 

Vfc  - Vph 

where  v*  is  the  average  velocity  of  the  beam,  v,/,  ^ is  the  root-mean-square  (rms)  spread 
in  the  axial  velocity,  and  Vpi,  is  the  phase  velocity  of  the  fastest  growing  ponderomotive 
wave.  When  the  axial  velocity  spread  is  large  enough  so  that  5 ;:g>  1,  the  FEL  operates 
in  a kinetic  regime  where  the  cold  beam  approximation  is  invalid.  This  generally  results 
in  a drastic  reduction  in  the  efficiency  of  the  FEL  since  only  a fraction  of  the  electrons 
contribute  to  the  growth  of  the  radiation. 

The  axial  velocity  spread  on  the  beam  is  initially  determined  by  the  beam  emittance 
and  energy  spread.  As  the  beam  propagates,  both  the  wiggler  fields  and  the  beam’s 
own  self-fields  cause  the  axial  velocity  spread  to  evolve.  When  the  beam  current  is  low, 
the  self-fields  can  be  neglected,  and  the  axial  velocity  spread  only  evolves  because  of 
the  betatron  oscillations  caused  by  the  wiggler  field.  This  leads  to  an  axial  velocity 
spread  in  the  form  of  a shear  Here,  is  the  axial  velocity  of  the  electron  fluid 
and  Vx  is  the  transverse  gradient.  When  the  self-fields  are  included,  they  oppose  the 
betatron  oscillations  and  reduce  the  shear.  Indeed,  in  the  case  of  a continuous  beam, 
an  equilibrium  can  be  found  such  that  forces  due  to  the  self-fields  exactly  cancel  the 
forces  due  to  the  focusing  fields  of  the  wiggler  [1].  In  this  paper,  we  perform  numerical 
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calculations  illustrating  this  equilibrium.  The  case  of  a pulsed  beam  will  be  considered 
in  a future  paper. 


NUMERICAL  MODEL 

We  study  the  electron  beam  flow  in  a PEL  by  numerically  following  the  orbits  of  a large 
number  of  particles  under  the  Influence  of  a prescribed  wiggler  field  and  a dynamically 
computed  self-field. 

The  particles  are  initially  loaded  into  a Gaussian  weighted  hyper-ellipsoid  in  six 
dimensional  phase  space.  They  are  advanced  in  time  using  the  standard  Boris  pusher 
[2]. 

The  self  field  model  is  derived  from  a method  used  in  the  Los  Alamos  code  TraceSD. 
First,  an  ellipsoid  is  fitted  to  the  particle  positions  via 

where  N is  the  number  of  particles,  i is  the  particle  index,  a is  one  of  the  three  spatial 
coordinates,  and  are  the  coordinates  of  the  beam  centroid.  The  force  exerted  on  a 
particle  by  the  self  fields  is  then  approximated  by 


= 5-3/2  ^ 30  1 / ^ 

4ti£o 

(1) 

^5-3/2  ? 30  1-/ 

4%eo'f  OyiOx  + Oy)C,, 

(2) 

f = 5-3/2^3qP_J_2 
47reo  OxOyOz 

(3) 

where  Eq  is  the  permittivity  of  free  space,  q is  the  charge  of  the  particle,  Q is  the  total 
charge  in  the  beam,  y is  the  relativistic  factor  of  the  beam,  and  / is  the  “space  charge 
form  factor”  which  was  tabulated  by  Lapostolle  using  data  obtained  numerically.  In  this 
paper,  ^ which  leads  to  / -^  0.  The  factor  of  5“^/^  converts  the  dimensions  of 
the  Gaussian  weighted  ellipsoid  into  the  dimensions  of  the  equivalent  uniform  weighted 
ellipsoid. 

The  wiggler  field  is  modeled  according  to  the  formula  given  in  Ref.  [3]  for  parabolic 
pole  pieces: 

B = Bq  sinh(2“'/^^ivA:)sinh(2”’'^^I:,t.y)cos(/:vvz)x 
-f  cosh  (2“'  ’ ^^kwx)  cosh  (2“  ’ ^^kwy)  cos{kwz)y 

cosh(2“ ' l^k^x)  sinh(2“ ' ^^k„y)  sin  {kwz)z 

Here,  k^  is  the  wiggler  wavenumber  and  x,y,z  are  the  unit  vectors.  The  electron  beam 
propagates  in  the  z-direction,  with  the  wiggle  motion  in  thex-z  plane.  The  parabolic  pole 
pieces  provide  equal  focusing  in  both  transverse  directions. 
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FIGURE  1.  Variation  of  spot  size  with  propagation  distance  for  a beam  with  and  without  space  charge. 
Solid  circles  represent  data  from  a calculation  with  J = 0.  open  squares  represent  data  from  a calculation 
with  i = 8 kA/cm^.  As  expected,  self  fields  extend  the  betatron  wavelength. 


SELF  COMPENSATION 

In  the  case  of  a continuous  beam,  the  effect  of  self  fields  is  to  modify  the  betatron 
wavenumber  according  to  [1] 

% = ^Po(I-SFP)'/2  (4^ 

where 

VYzoPzoYo  / 

Here,  = Anne^jmc^,  n is  the  electron  density,  e is  the  electronic  charge, 
^pQ  = a„,/t„,/\/2YoP;o  is  the  betatron  wavenumber  due  to  the  wiggler  field  only, 
= eBo/k^mc^  is  the  peak  normalized  vector  potential  of  the  wiggler  field,  p,-o  is 
the  axial  velocity  of  the  beam  normalized  to  c,  and  = (1  - Although  this 

formula  was  derived  for  flat  pole  pieces,  it  remains  valid  for  parabolic  pole  pieces 
as  well.  Note  that  if  the  self  field  parameter  (SFP)  is  unity,  the  betatron  wavenumber 
vanishes  and  the  beam  is  in  equilibrium. 

In  the  case  of  finite  emittance,  the  matched  beam  condition  becomes 

e«  = YoPzO^pa? 
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FIGURE  2.  Axial  velocity  distribution  for  a continuous  beam  10  cm  into  the  wiggler  field.  The  solid 
line  is  the  data  for  d = 0,  the  dashed  line  is  the  data  forj  = 8 kA/cm^ 


where  e„  is  the  normalized  emittance  (in  units  of  length  • radians),  and  Or  is  the  radius  of 
a uniform  density  continuous  beam.  For  a finite  spot  size,  equilibrium  occurs  for  > 0. 
However,  for  emittances  typical  of  a high  quality  photo-injector  [4, 5],  the  beam  is  nearly 
in  equilibrium  for  fcp  = 0. 

We  now  consider  numerical  calculations  of  the  electron  orbits  in  a wiggler  field  with 
Bo  = 3 kG  and  = In/k^  = 5 cm.  The  beam  had  initial  dimensions  of  0;t  = ctj,  = 1 mm, 
and  initial  emittance  of  zero.  Fig.  1 shows  measurements  of  Cy,  at  several  z positions.  The 
solid  circles  correspond  to  the  case  where  the  current  density,  J,  vanishes.  In  this  case 
the  beam  undergoes  betatron  oscillations  at  the  expected  frequency.  The  open  squares 
correspond  the  case  where  7 = 8 kA/cm^  which  gives  SFP  = 0.8.  In  this  case  the  spot  size 
is  nearly  constant.  According  to  Eq.  (4),  the  best  equilibrium  should  occur  for  SFP  = 1 
. In  the  numerical  calculation,  however,  SFP  = 1 led  to  a larger  deviation  in  spot  size 
than  SFP  = 0.8.  The  discrepancy  may  be  related  to  the  conversion  between  Gaussian 
and  uniform  beams. 

Fig.  2 illustrates  the  effect  of  space  charge  on  the  axial  velocity  distribution.  The  solid 
curve  corresponds  to  the  case  where  7 = 0 while  the  dashed  curve  corresponds  to  the 
case  where  7 = 8 kA/cm^.  The  distribution  is  significantly  narrower  in  the  case  of  the 
high  current  beam.  This  is  a result  of  the  fact  that  in  the  presence  of  strong  self-fields 
forward  momentum  is  not  lost  to  betatron  motions. 
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CONCLUSIONS 


Numerical  calculations  support  the  assertion  that  self-fields  can  improve  the  quality  of 
an  electron  beam  for  FEL  applications.  By  opposing  the  focusing  forces  of  the  wiggler 
field,  self-fields  can  increase  the  betatron  wavelength  and  reduce  the  spread  in  axial 
velocity.  This  reduces  S and  improves  the  performance  of  the  FEL. 
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